Uphill solitary waves in granular flows 
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We have experimentally observed a new phenomenon in the surface flow of a granular material. A 
heap is constructed by injecting sand between two vertical glass plates separated by a distance much 
larger than the average grain size, with an open boundary. As the heap reaches the open boundary, 
"soliton-like" fluctuations appear on the flowing layer, and move "up the hill" (i.e., against the 
direction of the flow) . We explain the phenomenon in the context of stop-and-go traffic models, and 
show that soliton-like behavior is allowed within a Saint- Venant description for the granular flow. 

PACS numbers: 45.70.-n,45.70Mg,45.70.Vn,81.05Rm,89.7S -W 



The rich dynamics of granular matter -studied for cen- 
turies by engineers- has attracted much attention from 
the Physics community since the early 1990's 000. 
Granular flows, for example, have concentrated intense 
interest, due to their relevance to natural avalanches and 
industrial processes, and also because they make liquid- 
like and solid- like behaviors coexist. They have been the- 
oretically described based on the existence of two phases: 
the rolling (or flowing) one, and the static one. Such 
idea has been casted into ad hoc phenomenological equa- 
tions 0, 0, @i mT, ° a Saint- Venant hydrodynamic ap- 
proach conveniently modified to take into account the 
particularities of granular matter 0, , and eventually 
by defining an orderparameter characterizing the local 
state of the system [fj, |lTjj . Finally, granular flows have 
been described by "microscopic" equations based on the 
newtonian motion of individual grains submitted to grav- 
ity, shocks, and trapping events [llLll^. 

In particular, researchers have studied both exper- 
imentally and theoretically granular flows on inclined 
planes and tubes, and granular heaps, finding a whole 
jungle of patterns, spatio-temporal structures, and other 
nontrivial phenomena, such as fingerin g |l3l. avalanches 
extending both downhill and uphill |14l . " logitudinal 
vortices " flfjl Hr^ |. "bubbling flows" Il7l llSj. "revolving 
rivers" [ljj and even "singing" dunes |2fJ. 

So it looks very unlikely that one can find further un- 
explored phenomena in experiments as simple as pouring 
sand on a heap with an open boundary. However, we re- 
port here the existence of soliton-like waves moving uphill 
in such experiments, i.e., the spontaneous appearance of 
bump- like instabilities in the flowing layer that propagate 
uphill, contrary to the flow of sand. 

We used sand with a high content of silicon oxide and 
an average grain size of 100/im from Santa Teresa (Pinar 
del Rio, Cuba) [lij . The sand was poured into a cell 
consisting in a horizontal base and a vertical wall, sand- 
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FIG. 1: Simplified diagram of the experimental setup 

wiched between two square glass plates with inner sur- 
faces separated by a distance w in the range from 0.3cm 
to 3cm (Fig. 1). The lengths of the base and the ver- 
tical wall were approximately d ss 23cm. The sand was 
poured vertically into the cell near the vertical wall us- 
ing tiny funnels with several hole diameters, in order to 
obtain different flux values. As the sand was poured into 
the cell, a heap grew until it reached the open bound- 
ary, where the grains were allowed to fall freely (Fig. 
1). Digital videos were taken using a High Speed Video 
Camera Photron FASTCAM Ultima- APX model 120K 
in the range from 50fps to 4000/ps, with a resolution of 
1024 x 1024 pixels. 

Fig. 2 contain images from experiments with w = 
0.7cm and F = 0.6cm 3 /s. Fig. 2 (a) is a picture taken 
from a video recorded at 50/ps, and show a closeup of a 
section of the surface at the central region of the heap, 
where our basic findings are easily identified: "bump- 
shaped" instabilities appearing at random places of the 
surface move "up the hill" , as indicated by the white 
arrows. These bumps could travel either alone or in tan- 
dem. The bumps maintain their "identity" within vari- 
able times, which can reach more than 10 seconds. After 




FIG. 2: Basic experimental findings, (a) Uphill bumps moving in tandem, identified by arrows, (b) Difference between 
two pictures of a single bump, separated by 50ms. In both pictures, the smallest ticks in the scales correspond to 1mm. 
(c) Spatial-temporal diagram of bump dynamics, where x grows to the right and t grows upward, spanning 20cm and 2s, 
respectively. 



that, they flatten out. We have also observed that, if 
two of them move at speeds different enough to interfere 
with each other, they can be still identified after the in- 
teraction. However, these are rare events, and we have 
no video record of them. All in all, our observations hint 
at a soliton-like behavior in the observed instabilities. 

Fig. 2(b) shows the difference between two pictures 
from a single bump taken from a video recorded at 
500/ps, which are separated 50ms from each other. The 
darker horizontal band allows to visualize the downhill 
motion. Careful examination of picture 2(b) allows to 
separate an approximately lmm-depth band of "flying" 
grains, and an approximately 1mm depth layer of "flow- 
ing" grains underneath. The height of the perturbation 
above the unperturbed stream is difficult to measure, but 
it can be roughly estimated as 1/5 of the flowing layer. 

Fig. 2(c) shows a spatial-temporal diagram of the free 
surface. A horizontal line of the video record was taken 
just at the (average) free surface, so lower parts appear 
black, and higher parts appear white. From the picture 
it can clearly be seen four shock-waves, passing suddenly 
from black to white, moving upward with a well defined 
velocity (even if it fluctuates a little on the right). On 
the left the appearance of a new soliton can be seen, first 
increasing in height before starting suddenly to propa- 
gate. 

The average speed of the uphill motion of the bumps 
as a function of the input flux, F was estimated from 
the analysis of the motion of several instabilities using 
pictures analogous to those used to obtain Fig 2(b). Fig. 
3 presents the dependence of the average uphill speed 
versus input flux. The solid line corresponds to v up ~ 
F 1 / 4 , which will be discussed below. 

Fig. 4 shows a sequence of pictures of a bump sepa- 
rated by intervals of 0.125s, extracted from a video taken 
at 4000 fps. Careful inspection of the video associated to 
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FIG. 3: Uphill speeds of bumps measured as a function of 
input flux, for w = 0.7mm 



this sequence of images reveals the mechanism of move- 
ment of the instability at the single grain level. The 
grains involved in the flowing layer in picture (a) show 
two types of behavior: while most of the grains keep flow- 
ing from left to right at an average speed Vfi ow ~ 10cm/ s 
as the bump passes underneath from right to left , a frac- 
tion of the grains at the lower part of the flowing layer 
show a different behavior. The white circle in Fig. 4(a), 
represents one of those grains, which is moving from left 
to right. From 4(b) to 4(d), the grain has incorporated 
into the static layer (deposition), while the bump is pass- 
ing from right to left above it. In Fig. 4(e), the grain has 
reincorporated into the flowing layer (erosion), moving 
again from left to right. 

Exchange of particles between the static and fluid 
phases has been considered in several models of surface 
flows since the 1990's Q. We believe that our bumps 
form when a fluctuation implies an extra deposition of 
grains from the flowing layer on the static one, which then 




FIG. 4: Sequence of pictures of an bump moving uphill, 
with the same parameters described for Fig. 2(e). The white 
circle represents a grain of sand (a) moving downhill as part 
of the flowing layer, (b)-(d) trapped in the static layer as the 
bump passes by from left to right, and (e) moving again as 
part of the flowing layer. The horizontal length of the picture 
corresponds to 22mm in reality 

propagates backwards through an "stop-arid-go" mecha- 
nism observed in traffic dynamics |2^ . However, in our 
case the usual traffic models should be modified to take 
into account the free surface. A way to do it could be a 
two-lane system consisting in a fast lane where cars flow 
steadily, and a slow lane where cars stop and go, produc- 
ing a backward wave (or train of wa ves) within the slow 
lane, well described in the literature |22j ]. 

All in all, the situation can be pictured as a solitary 
wave crawling uphill underneath a shallow stream of sand 
flowing downhill. In fact, the stop-and-go traffic model 
has been mapped into the "classical" Kortweg-de Vries 
(KdV) equation after a few approximations, resulting in 
soliton-like solutions (2^. In the context of the KdV 
equation, the speed of the soliton can be estimated as 
v s ~ \J gho where g and ho are the acceleration of gravity 
and the depth of the unperturbed stream, respectively 
[2^|. If one assumes the well accepted result that the 
depth of the flowing layer is proportional to F 1 / 2 [2l| . 
we then get v s ~ F 1 / 4 , which follows the experimental 
result shown in Fig. 3. 

Knowing that the KdV equation can be derived from 
mass and momentum conservation equations applied to 
an incompressible fluid flowing down a shallow channel 
[2^|. we could expect that soliton solutions can be found 
if we describe our system through hydrodynamics Saint - 
Venant equations modified to take into account the par- 
ticularities of granular matter 0, @ which we are writing 
relative to a reference frame parallel to the average free 



d t ( + d x (TH 2 /2) = (1) 

where £ is the height of the free surface, H is the thick- 
ness of the flowing layer, and V is the velocity gradient 
transversal to the flowing layer, which is assumed con- 
stant. Then, the momentum q = TH 2 /2 evolves with H 
due to the erosion/deposition process, like the one illus- 
trated for one grain in Fig. 4. The second equation cor- 
responds to the conservation of momentum, which can 
be seen as an equation governing the evolution of the 
thickness of the flowing layer, H: 

d t H + d x {TH 2 /2) = ^(tan9 - f i(H)) (2) 

where (J-(H) is the friction acting on the layer and 
tanO — —d x ( is the slope of the free surface. 

Let us assume that that the friction coefficient depends 
on H in the general form: 

M = a + 09 + ^)(h 2 - H 2 )^5~e{h 2 -H*) (3) 

where the values of constants a, /3, 7, S and e can be 
expressed in terms of experimental parameters v up , Vfi ow , 
ho, g, and a characteristic length Lq. Then, equations 
(1) and (2) can be transformed into the "textbook" KdV 
equation [24J 

LoVfiovTi" + 6(v up - Vfi ow )r) + 30 Vf l° w r) 2 = (4) 

no 
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In (4), 77 = C - h and 77" = where <p = x + v up t. 
One particular solution of that equation, is 

C ^ ho + ^ sech 2 { ^El ) (5) 
5 Lq 
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where we have taken Lq = 10ho, Vfi ow — iTho = 
y/gho, and v up = \vfi ow in order to match our experi- 
mental observations (i.e., mainly that the uphill speed is 
roughly 1/3 of the flow speed, and the height of the per- 
turbation is approximately 1/5 of the thickness of the un- 
perturbed flowing layer). Formula (5) describes a "text- 
book", bell-shaped soliton moving upward. This solution 
holds for a friction coefficient given by equation (3) which 
is a positive, decreasing function of H when reasonable 
experimental parameters are introduced, in agreement 
with previous models 0. Although we do not discard 
that "non-textbook" solitonic equations might describe 
the observed bumps without imposing such strong con- 
straints, we stress that our objective here is to show that 
"textbook" solitons can be derived if the system displays 
an appropriate n(H). 

We now present a number of additional experimen- 
tal facts relative to our uphill soliton waves. First, they 
have been observed in the whole range of experimen- 
tal parameters explored, i.e. 0.15cm 3 /s < F < 3cm 3 /s 
and 0.3cm < w < 3cm. As many phenomena in gran- 
ular matter, the uphill waves depend, to some extent, 
on the way the heap is prepared, and show some degree 
of "memory" . When the heap is prepared from scratch, 
uphill bumps nucleate at random places near the center 
of the heap, and, within a few seconds, they typically 
appear near the open boundary. If the experiment is 
stopped and re-started in such conditions, solitons reap- 
pear immediately. That happens even if a layer of sand 
less than lmm-thick is removed from the surface. If a 
thicker layer is removed, one has to wait a few seconds to 
observe them. This is probably due to the formation of a 
"compactified" layer during the flow that is necessary for 
the formation of these "jammed" bumps. Finally, all the 
results we have presented here have been observed just 
for one type of sand that produces sandpiles through "re- 
volving rivers" when dropped on a flat horizontal surface 
[lij. This fact supports the idea that a quite specific 
fJ,(H) dependence is needed to observe the solitary waves 
we report in this paper. 

In summary, we have observed soliton-like instabilities 
in a flow of sand established on a heap with open bound- 
aries, moving against the direction of the flow. The phe- 
nomenon can be understood in the light of stop-and-go 
traffic arguments, even for this free surface flow, and can 
be described by Saint- Venant equations adapted to gran- 
ular flows. The "microscopic" mechanisms that make a 
certain sand more suitable to show uphill bumps remains 
a mystery, but it seems clear that the dependence of the 
friction coefficient of our particular sand on the flowing 
depth is appropriate for the appearance of solitons. 
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